Singularities of Quantum Control Landscapes 

Rebing WuJ*| Jason Dominyij] Tak-San Hoj*] and Herschel RabitJ*| 

(Dated: July 14, 2009) 

Abstract 

A quantum control landscape is defined as the objective to be optimized as a function of the control 
variables. Existing empirical and theoretical studies reveal that most realistic quantum control land- 
scapes are generally devoid of false traps. However, the impact of singular controls has yet to be inves- 
tigated, which can arise due to a singularity on the mapping from the control to the final quantum state. 
We provide an explicit characterization of such controls that are strongly Hamiltonian-dependent and 
investigate their associated landscape geometry. Although in principle the singularities may correspond 
to local traps, we did not find any in numerical simulations. Also, as they occupy a small portion of 
the entire set of possible critical controls, their influence is expected to be much smaller than controls 
corresponding to the commonly located regular extremals. This observation supports the established 
ease of optimal searches to find high-quality controls in simulations and experiments. 



'Department of Automation, Tsinghua University, Beijing, 100084, P.R. China 

T Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA 
^Department of Chemistty, Princeton University, Princeton, New Jersey 08544, USA 

1 



I. INTRODUCTION 



Ma 



4|] was developed to evaluate the 



The concept of a quantum control landscape 
complexity of finding optimal controls, especially with respect to understanding the observed 
ease of achieving laser control of quantum systems flQ]. The search for optimal controls had 
been expected to be extremely difficult due to the complexity of quantum dynamics phenomena. 
However, practical studies show that such searches converge rapidly to high-quality solutions in 
optimal control theory (OCT) simulations as well as physically acceptable solutions to optimal 
control experiments (OCE) where many additional factors can be involved. 



A quantum control landscape [2, 



4|] is defined as the objective: 

7[6(-)] = fgkd) (l) 



at some given final time T, which is an implicit function of control field e(t) that steers the 
quantum system state if/(t) in an Af-dimensional Hilbert space H = C N through satisfaction of 
the Schrodinger equation: 

^(f) = [H + e(t)Him\ «K0) = <Ao, (2) 
at 

where the Planck's constant has been set to fi = 1. The free and control Hamiltonians, re- 
spectively H and Hi, are skew-Hermitian operators on the Hilbert space. The goal is to seek 
max e (.) /[£(•)], which entails a search on the landscape aiming to find at least one control that 
reaches the absolute global maximal value of J. 

The efficiency of searching for optimal controls is largely determined by the topology of the 
entire set of landscape critical points among which reside the ultimate desired optimal controls. 
Formally, a critical point corresponds to a control e(-) that satisfies 

6J=(VF( l /,(T)),6ifr(T)) = 0, V 5e{-), (3) 

where the inner product on C N is defined as (v, w) = ^(Vw). The corresponding state variation 
8i)j(T) at t = T is implicitly dependent on the control variation 8e{-) around e(-). 



In previous studies, we showed that when (i) the system is controllable at t = T, and (ii) any 
admissible control is regular (i.e., the mapping from de(-) to Si//(T) is surjective), the landscape 
is topologically equivalent to that of 

m = fw) (4) 

on the unit sphere of < H. Such landscapes are called kinematic in the sense that their topology 
is independent of the dynamics. The study of several classes of quantum control problems 
IS l3l |4t lZL [8|], including observable expectation- value optimization and quantum gate fidelity 
optimization, revealed that in principle no traps (i.e., local suboptima) exist to impede the search 
for optimal controls, thereby providing strong support for the observed ease of finding globally 
optimal controls in simulations. In the laboratory, optimization is also very efficient even for 
highly complex systems, but constraints on the controls likely imply that less than the absolute 
maximum value of the landscape is actually reached. 

In general, satisfaction of the controllability assumption is reasonable under generic circum- 
stances JsIIk]] because the Lie algebra rank condition for quantum controllability [9, 11 ] is easy 



to fulfill. However, the regularity of admissible controls calls for careful assessment because 
singular controls may exist corresponding to some VF(ip(T)) + without violating the condi- 
tion ©. In such a case, the criticality of a control is caused by its singularity. Understanding 
whether or not such controls are locally maximal (i.e., false traps) is important in obtaining a 
complete understanding of quantum control landscapes. 

In optimal control theory, there is evidence that singular controls may become local optima 
or even global optimal solutions (e.g., in time optimal control of rockets Jl^L . For quantum 



Hi 



systems, very few related studies seem to exist. Boscain and Chariot [13] proved that singular 
controls cannot be critical for a class of quantum systems that have multiple independent control 
interactions under the rotating wave approximation ; D'Alessandro showed that at most one 
singular control can be critical for the minimal-fluence control of two-level systems and the 



control is always constant II 1411 . For multi-level systems, there appears to be no general results. 



This paper will give a characterization of singular controls of single-input quantum systems 



and investigate their impact on quantum control landscapes. The balance of this paper is ar- 
ranged as follows. Section II defines the singular controls and Section III provides methods to 
compute singular controls. In Section IV, the landscape critical points are classified and com- 
puted in numerical simulations, and their impacts to the landscape are discussed in Section V. 
Finally, conclusions are presented in Section VI. 



II. CHARACTERIZATION OF SINGULAR CONTROLS 



From the viewpoint of functional analysis Ill5h . a control is singular (resp., regular) if the 
Frechet derivative 

d£^°' r : Se(-) h> 6iff(T), (5) 

of the end-point mapping E^ 7 : e(-) i-> iff(T) (defined on a neighborhood of e(-) for the 
L°°[0, T] norm) is rank deficient (resp., surjective) from the tangent space of e(-) to the tangent 
space T^t) of the unit sphere of *K at if/(T). The corresponding i^(-) is called a singular (resp., 



regular) trajectory of ©■ It may be shown lllql that E^ 0,T is Frechet differentiable with respect 
to the 1} topology on [0, T] (and therefore also with respect to the L°° topology. As a conse- 
quence, we can perturb © along the reference trajectory (driven by the reference control e(-)) 
to obtain an explicit form of the derivative: 

j t m) + sm) = {h + [6(o + 6e(t)m}m) + m)\ 

and then reduce it into the following time-dependent linear system after omitting higher-order 
terms 

^&K0 = A(t)5ifj(t) + B(t)6e(t), 6i/r(0) = 0, (6) 

where A(t) = H Q + e{t)H\ and B(t) = Hi^/(t). Let U(t) be the system propagator that evolves 
i//(0) to i]/{t) = U{t)i]/{0), then integrating © gives an expression for the derivative: 

6if/(T) = U(T) f U\t)H x ilj{t)56{t)At = U{T) f H x {t)i]/ Q 8e{t)At, (7) 
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where Hi(t) = W{t)HiU(t). Let M be the number of linearly independent functions of time 
among the real and imaginary vector components of Hi(t)if/ . By the above expression, a control 
is regular if there are M = IN - 1 (i.e., the dimension of unit sphere in *H) linearly independent 
functions over [0, T]. Otherwise, the control is singular and the number k = 2N - 1 - M is 
called its corank. 

In quantum optimal control theory, the cost function often appears in terms of the system 
propagator, i.e., J(e(-)) = F(U(T)) (e.g., F(U) = \Tr(W*U)\ for maximizing the gate fidelity 
with some specified unitary W), where U(t) obeys the evolution equation 

%-U(t) = [H + e(mW(t), U(0) = I N . (8) 
at 

For such problems, the singular controls correspond to those such that the Frechet derivative 

dE* oJ : Se(-) i-> SU(T), (9) 

of the end-point mapping E^ hT : e(-) i-» U(T) (defined on a neighborhood of e(-) for the 
L°°[0, T] norm) is rank deficient from the tangent space of e(-) to the tangent space 7~u(t) of 
the unitary group tl{N) at U(T). One can similarly derive the Frechet derivative from 8e(-) to 
SU(T) at U(T) 

6U(T) = U(T) \ H x (f)8e(t)&t, (10) 
Jo 

according to which a control is singular to the control-to-propagator mapping if H\(f) contains 
less than Af 2 linearly independent functions. For example: 

Example 1 Any constant control is singular to the control-to-propagator mapping. Their 
coranks are at least N - 1. 

Proof: Let e(t) = c be a constant control, then its singularity is equivalent to the linear depen- 
dence of the matrix elements of 

Hi(t) = exp(-t(H + cH{))Hi exp(t(H + cH{)) 
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as functions of time. Suppose that H + cH\ = QAQ^ where A is diagonal and Q is a unitary 
transformation, then the analysis is equivalent to investigating the matrix elements of 

Si(t) = &H x (f)Q = expHA)^ exp(?A), H 1 = $H X Q. 

The N diagonal matrix elements of H\{t) are all constant, implying that they are mutually 
linearly dependent, thereby the corank is at least N - 1 . The corank increases when Hi has 
fewer nonzero eigenvalues or H + cH\ has a degenerate spectrum. Q.E.D. 

It should be noted that the condition for a control to be singular to the control-to- state map- 
ping is stronger than that to the control-to-propagator mapping, i.e., any singular control (e.g., 
the constant control) for the control-to- state mapping must also be singular for the control- 
to-propagator mapping, as can be easily seen from © and (flOl) , but the inverse is not true. 
There exist singular controls for the control-to-propagator mapping that are not singular for the 
control-to-state mapping. 



HI. COMPUTATION OF SINGULAR CONTROLS 

This section will provide two approaches to numerically compute singular controls from 
different perspectives. The first one derives the singular controls by projecting the singular 
trajectories from a lifted space. The second one directly gives the control in a feedback form, 
however, additional smoothness constraints on the control are posed. 

Firstly, according to © and ©, a control e(-) is singular if and only if there exists a nonzero 
vector (p T e T^j) sucn that 

(<p T ,U(T)U^t)H l m) = (uw\T)<p T ,H l m) = ^ V*e[0,r]. (11) 
By defining the conjugate vector <p{t) = U(t)W(T)<p T e 7~^ w , then we have 

lifr(t) = [H Q + e{t)Hmt\ <A(0) = <Ao, (12) 
at 

%<P(t) = [H + e(O#i]0(O, <KT) = fa, (13) 
at 
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subject to the algebraic constraint from (TTTI) 



(0(0,^1^(0) = o, Vfe[0,r]. 



(14) 



This is a two-point boundary-value problem in time which has to be solved by iterative numer- 
ical algorithms. As will be seen later, these equations can be also derived from the Pontryagin 
Maximum Principle with respect to a given cost function, with 4>(T) assigned to be the gradient 
vector of the cost function. 

Using (fT2l) and (PT3T) . the equation (fl4l) can be differentiated to derive a new algebraic con- 
straint: 



which can be again differentiated to arrive at an explicit relationship between a singular control 
f(-) and the corresponding state trajectory: 



—(fctXH^it)) = => (0(0, [H , [//o,#i]]<KO> + e(t)(<f>(t), [H u [# ,#i]]<K0> = 0. (16) 



According to (fT6l ), we classify singular controls as follows: 

(i) When (4>(t), [H\, [Hi, Ho]]tf/(t)) ± 0, Vf e [0, T], the singular control can be expressed in 
a feedback form 



One can then combine (PT2l) . (PT3T) and (PT71) to solve for the singular control. Such controls 
are called minimal-order singular controls. Notice that since each nonzero 4>t must uniquely 
correspond to some 0(0) = 0o ^ at ? = 0, we can equivalently integrate the differential 
equations (PT2l) and (PT3T) from f = (this can be done from the other end t = T as well) for any 
given pair of (i]/q, 0o) that satisfy 



Moreover, since (PT2l) and (TT3T ) share the same evolution propagator U(t), we can use the dy- 
namics of the system propagator: 



-Mt),H^(t)) = => <0(O,[#o,#i]<KO> = O, 



(15) 



e(t) = - 



(0(0, [gpJg^HijjjKO) 

[#i,[# ,#i]¥(o> 



(17) 



(<f>o,Hnf/o) = <0o, [#o,#i¥o> = 0, <0o, [#i, [#o,#i]¥o> * 0. 




(18) 
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to obtain a singular trajectory in the unitary group U(N) parameterized by (if/ , </> ). In this way, 
singular controls can be systematically generated without iterative computations. 

(ii) When (0(0, [H lt [HuHoWit)) = on [0, T], this implies (cf>(t), [H , [H u # ]]^(0> = 
as well. In such cases, equation (fT71) is not sufficient for determining a singular control. How- 
ever, one may go on differentiating these two quantities until e(t) can be explicitly expressed. 
Let H av .. ai = [H ai , [H a2 , [■■■ , [H ak { , H ak \ ■ • • ]]]. If there exists a finite integer k > 2 such that 
{<p{t),Hp v ~i3 { ty{t)) - for an Y C^i, • - • ,Pt) e {0, if, 2 < I < k, but the differentiation of some 
<0(f), H au ^(t)) gives 

<<K0, H 0av .. ak m) + eMWl H lav .. ak ifr(t)) = 0, (19) 

where (0(f), H laiv .. ak ifr(i)) ± 0, then a singular control can be formally obtained from a feedback 
equation 



At \ (fo.WW, U(t)iff ) 
and we call it a k-th order singular control. 

Let S m = spanj/Zi} and S (f) = span{H av .. a( \ai, . . . , cti = 0, 1} (£ > 2) be the subspaces of 

skew-Hermitian matrices generated by i tuples of commutations. A geometrical interpretation 

for a control to be fc-fh order singular is that the adjoint vector (p(t) e T^t) is orthogonal to 

B (e) i]/(t) c T m for any 1 < t < k but not orthogonal to B (k+l) i]/(t), i.e., (p(t) belongs to the 

following kth order singular cone at if/(t): 

'V® =Le T m I 1^, [J S (f V(o| = 0, (0, S ( * +1 ty(0> * j . 

To locate a fc-th order singular control, one can choose a pair {\f/ Q , <p Q ) such that O 6 < T> > an d 
integrate (1201) from t = 0, provided the solution exists and is unique. 

(iii) When the control function does not explicitly appear in (fT9~l) for any integer ieN, then 
its order is infinite. Let X be the Lie algebra generated by H and Hi, and Xo = U7=i ^ {t) 
be the minimal ideal in X that contains H\. The codimension of Xo<Ao m -£ i s either or 1. 
For infinite-order singular controls, the codimension of Xo^o i n -£<Ao must be 1 (otherwise <p(t) 

8 



has to vanish), i.e., the adjoint vector (p(t) e varies in the one-dimension complimentary 
subspace of -£oiA(0 m X^(0- 

It is also possible for the denominator in (fT71) to cross zero at isolated time instants, which 
divide a singular control into pieces of singular "arcs" whose orders may be different with each 
other. In this paper, we only consider singular controls whose order is constant on [0, T]. 

An alternative approach to produce (i) and (ii) above is as follows. Denote the 2N - 1 in- 
dependent elements in the vector Hi(t)if/ by = (£\(t), ■ • ■ ,£ 2 at-i(0)- For a singular control, 
there must exist a nonzero constant vector O = (A> ' 1 ' ,C2N-i) such that Ym=\ 1 c &(f) — 0, 
Vfe [0, T], which can be repeatedly differentiated to give 

2N-\ 

YjCi^W^O, k= ,2N-2. 



So we have 



<52N-\"> 



C2 



= 0, Vf 6 [0, T] 



) \ C2N-\ ) 

which implies that the Wronskian must vanish, i.e., 



&lV-l(0 



C-f(o 



= o. 



From the previous derivations, £ ; - 2) (0 can be decomposed into two parts as £;f\t) = 



af\t) + e(t)bY'(t), where a t (t) and bi(t) correspond to linearly independent functions in 



,(2) 
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[H , [H , Hi]](f)^o and [H U [H , //i]](0<Ao, respectively 



af(t) 



af(t) 



£ {2N - 2) (t) 

S2JV-1 W 



6(0 
^(0 
(0 



z> (2) 



6(0 



t2) (Z) 



(1) (0 



f (2W-2) () 



= 0, (21) 



where the time derivatives of e(0 (up to the (2iV - 4)-th order) are involved in the fourth to 
the (2N - 2)-th rows. Hence (I2TI) forms an ordinary differential equation (up to (2N - 4)-th 
order) of e(0 whose coefficients are functions of the system propagator U(t). This relation can 
be computed in the Schrodinger equation to solve for the singular controls from given initial 
values of e(t) and its derivatives. 

Both the above two approaches calculate singular controls via ordinary differential equations 
from the same condition ©. The solution by the first approach is parameterized by a prescribed 
vector 0o as the initial condition of the conjugate equation, while the latter is by the initial 
conditions of the time derivatives of 6(0- Moreover, it is easy to see that the order of the 
differential equation (f2TT) with respect to e(t) is IN - 2 - k, as the derivatives of 6(0 start to 
appear from the (k + 2)-th row, where k is the order of the singular control defined in the first 
approach. In this regard, these two approaches are equivalent. In comparison, the first approach 
is numerically more efficient and will be adopted in the simulation examples below. The latter 
is conceptually simple because it appears directly as the sum of the derivatives of the control 
function and the state without introduction of any conjugate vector, and hence provides a useful 
perspective for the origin of singularity. 



IV. SINGULAR CONTROLS AS LANDSCAPE CRITICAL POINTS 

As analyzed above, the critical points for a given control landscape can be a regular or 
singular control. The corresponding kinematic gradient must vanish when the critical point is 
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a regular control , but it may not vanish when the critical point is a singular control. In this 
regard, we can classify the landscape critical points into the following three categories: 

Definition 1 A control is said to be regularly (singularly) kinematic if it is regular (singular) 
and VF(if/(T)) = 0. Otherwise, ifVF(if/(T)) ± 0, it is said to be non-kinematic. 

Let Cf = ty(T) 6 S<h | VF(if/(T)) = 0}, where S<h is the unit sphere in < H, be the set of 
kinematic (either regularly or singularly) critical points. Any control that steers the trajectory 
onto Cp at t = T must be a kinematic critical control. As indicated in previous studies [1], the 
dimension of regularly kinematic controls is infinite and its codimension in the set of admissible 
controls is the same as that of C F in < H. 

For nonkinematic critical points, their corresponding kinematic gradient vector <f> T = 
VF(i//(T)) must be nonzero and belong to some singular cone < VS r) (2 < k < oo). This criterion 
forms a nonlinear constraint on the final state if/(T) and thereby defines a subset of ( H: 



which will be called the &-th order singular surface for the control landscape F . Every £-th 
order singular control has to terminate at this surface to become a nonkinematic critical point. 
Starting from an arbitrary point if/ T e S ®, one can determine a constant &-th order nonkinematic 
critical point by integrating (fT2l) and (fT3l) backwards in time along any proper direction in 
if the solution exists. 

It is difficult to estimate the dimension of the entire set of nonkinematic controls, as they 
can be a combination of singular arcs with different orders, which correspond to an infinite 
number of possibilities. Nonetheless, as illustrated in Fig. \T\ the set of regularly kinematic 
critical controls is much richer than nonkinematic ones because they may cross the surface Cf 
along any direction and with any admissible flows. By contrast, the set of singular (kinematic 
and nonkinematic) critical points is much more limited. 




(22) 
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FIG. 1: Schematic illustration for the geometry of the set of extremal controls. Every extremal trajectory 
ends up at the singular surface Sp at t = T. The singular trajectories (bold lines) follow the flow defined 
by (TT8T ) and traverse the singular surface along some direction (dash lines) in the singular cone, among 
which those that traverse Cf are singularly kinematic, and the remainder are nonkinematic. Any other 
control that steers the trajectory to traverse Cf is regularly kinematic. 

V. THE ROLE OF SINGULARITIES IN QUANTUM CONTROL LANDSCAPES 

We have derived that all of the nonkinematic critical points are subject to the equations 
(fT2l) and (fT~3b with the algebraic constraint (4>{t),Hiij/{t)) = 0, where the boundary condition 
satisfies (p T = VF(ijj(T)) + 0. On the other hand, any (regularly or singularly) kinematic critical 
points automatically satisfy the same differential equations and the algebraic constraint because 
4> T = VF(ifr(T)) = 0. Therefore, in the language of the Pontryagin Maximum Principle II 1711 . 
we can unify the conditions for a control to be a landscape critical point as: 

dib on dd> on an 

77 = ^7' -£ = -T7> 1T =0 ' <A(°) = , Ao, <Pt = VF(«KD), (23) 
at d(p at dip oe 

where H[e(f),if/(t),(f>(i)] = (<p(t),H ifr(t)) + e(t)((f>(f),Hii/f(t)) is the pseudo-Hamiltonian func- 
tion. 
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In common optimal control problems, the dynamics are often taken into account through the 
cost function 

J[e(-)] = F(i(r(T))+ f L(^(t), e(t),t)dt, (24) 
Jo 

where L(-) is a function chosen to balance the dynamical performance issues. The standard 
Pontryagin maximum principle then corresponds to the following pseudo-Hamiltonian function 
to solve for critical points of J (conventionally called extremals): 

H[«A(0, 0(0, e(t)] = AL(ifr(t), 6(f), + (0(0, «(0> + 6(0(0(0, #i<K0>, <KT) = VF(<KT)), 

where A is a constant. Taking the optimization process as a dynamical game between the end- 
point cost (the first term) and the dynamical part (the second term), the extremals corresponding 
to A + result from the trade-off between the two costs, under which the system cannot attain 



perfect yields II 1811 (i.e., the highest yield in the control landscape (QQ))- Such extremals are 
conventionally called normal extremals, which are generally not critical points of (QQ)- Normal 
extremals are always regular because the necessary conditions (PT21) and (fT3l) are not satisfied. 

Extremals corresponding to A = are called abnormal, and they are also critical points of 
(OQ). Under such controls, the end-point cost completely overwhelms the dynamical cost so that 
the resulting trajectories are independent of the choice of the cost function L(-) in the integral. 
They include all kinematic and non-kinematic critical points discussed in this paper (there is no 
analog of the kinematic picture for normal extremals because the dynamics is always relevant). 
Thus, we can classify the extremal controls as in Fig. [21 




FIG. 2: Classification of extremal controls for general optimal control problems. 
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The optimality of abnormal extremals can be analyzed through the second-order variation 
of®: 

6 2 J = OkdC W), W)) + <V7((A(r)), 5 % ^{T)), (25) 

where the kinematic part Q^t) is the Hessian quadratic form at if/(T). The second term vanishes 
when the control is a kinematic extremal, either regular or singular, and hence leaves the op- 
timality determined solely by the positive-definiteness of the kinematic Hessian form Q^t) on 
the set of achievable Sifr(T) at if/(T). For regular abnormal extremals (i.e., regularly kinematic 
critical points), the optimality is exactly the same as that of the kinematic Hessian, and hence 
their topology is reflected by that of the corresponding final state iff(T) in the kinematic picture, 
exhibiting a universal Hamiltonian-independent property. 

If an abnormal extremal is a singularly kinematic critical point, then its Hessian form Q^t) 
is defined on the range of the Frechet derivative AE^ 1 as a proper subspace of T^j)- When 
the unrestricted Q^t) is positive (negative) semi-definite, the restricted Q^t) is also positive 
(negative) semi-definite. However, in the case that the unrestricted Q^t) nas both positive and 
negative eigenvalues, its restriction on the range of d£*° r may be positive (negative) when 
the range of d£^° T is contained in the subspace spanned by the eigenvectors corresponding 
to positive (negative) eigenvalues. This implies that a saddle point may degenerate to a local 
minimum (maximum). 

Nonkinematic critical points are beyond the scope of the kinematic picture, for which the 
second-order variation 8 2 ij/{T) related to the system dynamics is nonvanishing. Hence, the 
optimality of nonkinematic critical points is much more complex to assess as the Hessian form 
has to be discussed on an infinite dimensional space of control fields. As indicated by Bonnard 



and Chyba Hl5H . there is a possibility that such critical points are local optima in the control 
landscape. However, in our numerical simulations, no such traps have been found. Figj3]shows 
examples of nonkinematic extremal controls for the quantum state transition control landscape 
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J[e(')] = Wfft(T)\ 2 for a four-level quantum system, where 



-0.50 



0.00 



0.20 



0.60 



0.30 0.75-0.20/ 0.65 0.40 

0.75 + 0.20/ 0.70 0.70 - 0.50/ 0.20 + 0.30/ 

0.65 0.70 + 0.50/ 0.30 0.50 

0.40 0.20-0.30/ 0.50 0.60 



with the initial and target states being ifr(0) = [1 0] r and if/f = [0 00 l] r , respectively. 
Singular trajectories were found by solving equation ([TBI) . After locating such singular controls, 
we started a gradient search from small neighborhoods of the controls, and found that the 
search always climbed towards perfect yield (7 = 1.0) without being trapped. This behavior is 
consistent with the observation that singular controls have not been located when performing 
common optimal control simulations, i.e., a gradient flow trajectory is not attracted to a singular 
trajectory even when one is nearby. 



VI. CONCLUSION 

This paper considered the role of singular controls upon the search for optimal solutions over 
quantum control landscapes. We indicate that the regularity assumption of admissible controls 
can be violated, which gives rise to singularly kinematic or non-kinematic critical points that 
are beyond the scope of the kinematic picture. For single-input systems, these singular controls 
could possibly be locally optimal. However, such cases were not found in our simulations. 
Moreover, since the entire set of singular controls is small in contrast with that of the regular 
controls, the overall landscape is not expected to be rugged. Hence, regular controls dominate 
over the singular ones, and the regular controls should still determine the overall landscape 
complexity. This conclusion is fully consistent with the general quantum control literature 
where no specific evidence is seen for singular controls to be local sub-optimal solutions. 
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FIG. 3: Nonkinematic critical points of the quantum state transition control landscape J[e(-)] = 
\i/Ai/f(T)\ 2 for a four-level system. The figure on the left displays two singular extremal controls 
from randomly chosen initial states. Starting from some control ej n itiai(0 in a small neighborhood 
(IkinitiaiO) _ e(OII ^ 0.01) of the singular control So(f), the steepest ascent algorithm is performed to 
search for optimal controls. Both of them climb and are not attracted by the corresponding singular 
controls, but approach a perfect yield, showing that these singular controls are not local false traps. 

When the system has multiple control fields associated with independent operators H^. 



the definitions of extremals can be extended but will not be given here. In this case the con- 
dition for a control to be singular and further critical is much more stringent. In particular, 



(H , Hi, • • • ,H m ), showing that almost all such systems do not possess minimizing singular 
controls. Therefore, the impact of singular controls on the control landscape is smaller yet, or 
even disappears, for multi-control quantum systems. 




(26) 



Chitour et al\\9] proved a generic property for the system class represented by (m + l)-tuple 



16 



Acknowledgments 



The authors acknowledge support from the NSF. JD acknowledges support from the Pro- 
gram in Plasma Science and Technology (PPST) at Princeton. Useful discussions with Profes- 
sor Bernard Bonnard are appreciated. 



[1] T.-S. Ho and H. Rabitz, J. Photochemistry Photobiology A 180, 226 (2006). 

[2] H. Rabitz, M. Hsieh, and C. Rosenthal, Science 303, 1998 (2004). 

[3] H. Rabitz, M. Hsieh, and C. Rosenthal, Phys. Rev. A 72, 52337 (2005). 

[4] R. Wu, H. Rabitz, and M. Hsieh, J. Phys. A 41, 015006 (2008). 

[5] H. Rabitz, R. de Vivie-Riedle, M. Motzkus, and K. Kompa, Science 288, 824 (2000). 

[6] V. Bonacic-Koutechy and R. Mitric, Chem. Rev. 105, 1 1 (2005). 

[7] T.-S. Ho, J. Dominy, and H. Rabitz, Physical Review A (Atomic, Molecular, and Optical Physics) 
79, 013422 (pages 16) (2009), URL |http : //link, aps . org/abstract/PRA/v79/eQ13422| 

[8] R. Wu, R. Chakrabarti, and H. Rabitz, Physical Review A 77, 052303 (pages 13) (2008). 

[9] V. Ramakrishna, M. V. Salapaka, M. Dahleh, H. Rabitz, and A. Pierce, Phy. Rev. A 51, 960 (1995). 
[10] R.-B. Wu, T.-J. Tarn, and C.-W. Li, Phys. Rev. A 73, 012719 (2006). 
[11] W. Boothby, J. Diff. Equa. 17, 296 (1975). 

[12] S. Serovaiskii, Counterexamples in optimal control theory (Brill Academic Publishers, The Nether- 
lands, 2004). 

[13] U. Boscain and G. Chariot, ESAIM: Control, Optimisation and Calculus of Variations (COCV) 10, 
593 (2004). 

[14] D. D'Alessandro, IEEE. Trans, on Automatic Control 46, 866 (2001). 

[15] B. Bonnard and M. Chyba, Singular trajectories and their role in control theory, vol. 40 of Math- 

ematiques and Applications (Springer, Berlin, 2003). 
[16] J. Dominy and H. Rabitz, Journal of Physics A: Mathematical and Theoretical 41, 205305 (21pp) 

17 



(2008). 

[17] R. Stengel, Optimal control and estimation (Dover Publications, New York, 1994). 
[18] W. Zhu and H. Rabitz, The Journal of Chemical Physics 109, 385 (1998). 
[19] Y. Chitour, F. Jean, and E. Trelat, SIAM J. Control Optim. 47, 1078 (2008). 



18 



